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Abstract 

The data for the mean squared transverse momentum {p1)[Et) as function of transverse 
energy Et of J and ip' produced in Pb-Pb colhsions at the CERN-SPS are analyzed 
and it is claimed that they contain information about the time structure of anomalous 
suppression. A transport equation which describes transverse motion of J/ip and tp' 
in the absorptive medium is proposed and solved for a QGP and a comover scenario 
of suppression. While the comover approach accounts for the data fairly well without 
adjusting any parameter, the fit to the data within the QGP scenario requires to assume 
anomalous suppression to become effective rather late, 3 — 4 fm/c after the nuclear overlap. 

The discovery in 1996 of anomalous J/ip suppression in Pb-Pb collisions at the SPS has 
been one of the highlights of the research with ultrarelativistic heavy ions at CERN . Does it 
point to the discovery of the predicted quark gluon plasma (QGP)? Six years later the situation 
is still confused, since several models - with and without the assumption of a QGP - describe 
the observed suppression, after at least one parameter is adjusted. The data on the mean 
squared transverse momentum {pf)[Et) j2] for the (this symbol stands for J/tp and tp') in the 
regime of anomalous suppression and as a function of transverse energy Et have received less 
attention - for no good reason. We claim: {p1){Et) contains additional information about the 
nature of anomalous suppression and may help to distinguish between different scenarios. In 
this paper we investigate how the time structure of anomalous suppression influences the values 
{p1){Et). This idea has already been considered more than 10 years ago 3J (c.f. also more recent 
works jHE]) and is based on the following observation (Fig^: Anomalous suppression is not an 
instantaneous process, but takes a certain time depending on the mechanism. During this time 
T/^'s produced with high transverse momenta may leak out of the parton/hadron plasma and 
escape suppression. As a consequence, low pt are absorbed preferentially and the {pf) of the 
surviving (observed) V^'s show an increase S{pf), which grows monotonically with the mean time 
tA, when anomalous absorption acts [Sj. In this letter we propose a general formalism of how to 
incorporate the effect of leakage into the various models, which have been proposed to describe 
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anomalous suppression and we extract information about the time tA from a comparison with 
experiment. 

It has become customary to distinguish between normal and anomalous values of suppression 
S{Et) = cr^ (Et) / a^^ (Et) for ip's produced in nuclear collisions, c.f. reviews [HIIZI- Here, 
and 0"^^ are the production cross section for a ip and a Drell-Yan pair in an AB collision, 
respectively. By definition, ip^s produced in pA collisions show normal suppression via inelastic 
ipN collisions in the final state and normal increase of (p^) (above {p1)nn in NN collisions) 
via gluon rescattering in the initial state. These normal effects are also present in nucleus- 
nucleus collisions and happen, while projectile and target nuclei overlap. Anomalous values of 
S and {pD are attributed to the action on the tp by the mostly baryon free phase of partons 
and/or hadrons (we call it parton/hadron plasma) which is formed after the nuclear overlap. 
It may lead to deconfinement of the ip via colour screening in the QGP, dissociation via gluon 
absorption or inelastic collisions by the comoving hadrons during the later period of the plasma 
evolution. In this letter we describe anomalous ip suppression within a transport theory and 
apply it to two rather different scenarios: (I) Absorption involving a threshold in the energy 
density like in a QGP scenario, (II) continuous absorption via comovers. 

We denote by da"^ /dpt{Pt, Et) the cross section for the production oi a. ip with given pt and 
in an event with fixed transverse energy Et. It can be related to the phase space density via 

-^{pt,Et)=\imJ dbP{Et;b) J ds f^{s,pt,t-b). (1) 

Here, P{Ef, h) describes the distribution of transverse energy Et in events with a given impact 
parameter h between projectile A and target B. We follow ref. |B] in notation for P{Et; b) and 
the values of the numerical constants. The function f^{s,pt,t]b) is the distribution of ip^s in 
the transverse phase space {s,pt) at time t for given b. 

We define t = as the time, when the process of normal suppression and normal generation 
of (pt) has ceased and denote by ft[{s,Pub) the distribution of ip's at this time. is taken 
as initial condition for the motion and absorption of the tp^s during the action of anomalous 
interactions. The evolution of the ip is described by a transport equation 

^^f + VfVsf = -af. (2) 

The time dependence arises from the free streaming of the ip with transverse velocity Vt = 
Pt/y^m^ + Pt (l.h.s.) and an absorptive term on the r.h.s., where the function a{s,pt,t;b) 
contains all details about the surrounding matter and the absorption process. We have left 
out effects from a mean field, because the elastic ipN cross section is very small and have also 
neglected a gain term on the r.h.s., because recombination processes c + c + N—>-tp + X seem 
unimportant at SPS energies where at most one cc is created per event. 
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Figure 1: Schematic picture of the leakage phenomenon. Between the Lorentz contracted rem- 
nants "B" and "A" of the nuclei which have collided charmonia move in the created par- 
ton/hadron plasma. Those ip's with large transverse velocities Vt (case (1)) may leak out and 
escape suppression, while low Vt particles may remain (case(2)) leading to an effective increase 
of {pD for the surviving (observed) ip's. 

Eq. Q can be solved analytically with the result 

=exp(^- J^dt'a{s-Vt{t~t'),pt,t';b)^fUs-Vtt,Pub), (3) 

which for t = reduces to f^ = f^. If we denote by t/ the time when anomalous suppression 
has ceased, a{s,pt,t; 6) = for t > tj, the limit t ^ cxd in eq. ((T)) can be replaced by setting 
t = tf, since the distribution in pt does not change for larger t's. 

There is little controversy about ip production and suppression in the normal phase: The 
gluons, which fuse to the cc, collide with nucleons before fusion and gain additional pf. The ip 
on its way out is suppressed by inelastic collisions without any change in (p'f). Neglecting 
effects of formation timel^l, one has 

ft{s,Pub) = aff^ J dzAdzB pA{s,ZA)pB(b- s,zb) ■ 

■ exp {-(^tbsiTAis, za, oo) + Tb(6 - s, -oo, zb)]) {Pt)^^ exp (-p^/ (p2^^^ ^ (4) 

where 

(j)t)N(b, s, Za, Zb) = {pI)%n + Sa^Po M^-a(s, -00, za) + Tb(6 - s, zb, +00)] (5) 
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with the thickness function T(s, ^1,^2) = Izi d-z p{s, z). All densities Pa, Pb are normalized to 
the number of nucleons (po is the nuclear matter density). We shortly explain eqs. (jH) and ©: 
For given values b and s*in the transverse plane the ip is produced at coordinates zj^ and zb in 
nuclei A and B, respectively. On its way out, the ip experiences the thickness Ta{s, za, 00) and 
TB{b — s, —00, Zb) in nuclei A and B, respectively and is suppressed with an effective absorption 
cross section (Tabs- The two gluons which fuse carry transverse momentum from two sources: 
(i) Intrinsic pt, because they had been confined to a nucleon. The intrinsic part is observable in 
NN — > ip coUisions and leads to {Pt)tiN ^q. (ii) Collisional contribution to pt, because 
in a nuclear collision, the gluons traverse thicknesses Ta{s, —00, za) and Tb(6 — s, zb, +00) of 
nuclear matter in A and B, respectively, and acquire additional transverse momentum via gN 
collisions. This is the origin of the second term in eq. (0). 

The constants cr'^lg , crths ^^^1 QgN are usually adjusted to the data from pA collisions, before 
one investigates anomalous suppression. Fig. 2 shows (dashed curves) the results for normal 
suppression S^{Et) and {p^t)^{Et) calculated with eq. © in eq. (0). While the difference 
between calculation and data is enormous for the suppression, it is rather small for {pf)'^{Et). 

Since the physical origin of anomalous suppression is not yet settled, we investigate sup- 
pression S^{Et) and {p1)^{Et) for two models, which have rather contradictory assumptions. 

I. Threshold (QGP-) scenario: -i/^'s are totally and rapidly destroyed, when they are in a 
medium with energy density above a critical one, and nothing happens elsewhere. As a 
representative model we use the approach by Blaizot et al. jH]. 

II. Comover scenario: The plasma of comoving partons and/or hadrons leads to a continuous 
absorption of long duration due to inelastic collisions with the comoving particles. As a 
representative model, we use the approaches by Capella et al. and Kharzeev et al. 

m 

In this letter we study the effect of leakage on the observed values of {pf) within two well 
established scenarios, using their assumptions and parameters. We do not introduce any mod- 
ifications like (i) the pt dependence of the absorption process (i.e. pt dependence of a{s,pt,t; h) 
in eq. (j21)) and (ii) expansion of the plasma during absorption. Both effects may contribute to 
{pD (we will present a qualitative discussion at the end), but both require a detailed study by 
themselves. 

We begin with model I: In their schematic approach Blaizot et al. include anomalous 
suppression via 

f^{sJt,h) = e{n,-n,{h,s))fUsJt,h). (6) 
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Here ric is a critical density and np(6, s) is the density of participant nucleons 

np{b, s) = Ta{s, -oo, +oo)[1 - exp [-a^fTeib- s, -oo, +00))] + {A ^ B). (7) 

According to eq. © all ip^s are destroyed if the energy density (which is directly proportional 
to the participant density) at the location s of the is larger than the critical density. All other 
V^'s survive. While the prescription eq.® successfully describes the data in the full Et range 
of anomalous suppression after the only one free parameter, n^, is adjusted, the predictions for 
{pf){Et) are significantly below the data, especially at large Et (see below). 

The expression eq. (jHl) for the phase space distribution including anomalous suppression 
within the threshold model is recovered within our transport approach eq. ^ by setting 

a{s,p,t; h) = ao9{np{b, s) — nc)S(t) (8) 

and taking the limit ao 00. The delta function 6{t) has to be included to recover the 
expression eq. © and may be understood by the physical picture that the energy density is 
highest at t = and anomalous absorption most effective then. 

There are various ways to introduce another time structure into the absorption term. We 
will try two options, one being 

a{s,pt,t;b) = ao9{np{b,s) - nc)S(t -Ia). (9) 

The idea of a threshold density is kept, but suppression does not act at t = but at a later 
time tA, which time is then determined from a comparison with the data. For times t > tA, 
one finds from the general solution eq. Q (and — > 00) 

f{s,pt,t)=e{n,-np{b,s))ff^is-VttA,Pt), t > Ia (10) 

which differs from the expression by the motion in phase space of during the time 
< t < tA- For t > Ia the momentum distribution derived from does not change any more. 

The suppression S'^{Et) and {pf)ip{Et) calculated with the distribution from eq. (fTn|l 
depends on the parameter tA- We use the parameters of ref. |H] where available, i.e. for J/V': 
^abf ~ 6-4 mb, ric = 3.75 fm^^. The parameters for the generation of (p^) by gluon rescattering 
are taken from a fit to the pA data [2 : {pDnn = l-H (GeV/c)^ and ag^ = 0.081 (GeV/c^) 
fm~^. We also account for the transverse energy fluctuations jH] which have been shown to be 
significant for the explanation of the sharp drop of J / ip suppression in the domain of very large 
Et values, by replacing Up by -j^np where {Et) is the mean transverse energy at given b. We 
then calculate a"^^"^ /a^^ as a function of tA- Since the critical density for J is quite high, the 
leakage affects only the very high momentum t/^'s. Since their number is small, the calculated 
results for the suppression S^/'f'{Et) depend only little on tA (Fig. 2). 
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We turn to a discussion of {p'l)'^/^{Et). Fig. 2 shows calculated curves for values of = 
to 4 fm/c. The dotted line {Ia = 0) is the result of the original threshold model with immediate 
anomalous suppression (and has been predicted in [IH])- It fails badly at large values of Ef. 
Also no other curve with a given Ia describes the data for all values of Et. We have to conclude 
that Ia depends on Et. tA{Et). The larger the values Et the later anomalous suppression 
acts. From a comparison with data we have ^ 2 fm/c for Et < 80 GeV, tA — 2.5 fm/c for 
80 <Et< 100 GeV, and ^ 3.5 fm/c for Et > 100 GeV. 
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Figure 2: Nuclear suppression /a^^ and {pf)^ for J/ip (above) and ip' (below) as a function 
of transverse energy Et. Data are from 0]/ for /a^^ and from ^ for (pf). Dashed curves 
show the result of normal suppression alone. The dotted lines correspond to the original thresh- 
old model, which in our notation corresponds to = 0. The other curves include anomalous 
suppression within the threshold model eqs.^ and I\1U\} . where anomalous suppression is as- 
sumed to act at time tA > 0. The curves are labeled by the values tA = 1,2,3,4/m/c for J/ip 
and tA = 2, 4, 6 fm/c for ip' . Also the curves in o"^ jo^^ carry these labels, lowest curve tA = 
and monotonic increase with tA. 



A similar analysis is performed for the tp'. The data for the suppression are taken from 
ITH] . those for {pf)'^' from Since the ijj' has not been treated by Blaizot et al., we fit 



the suppression data and find values for the absorption parameters, a 



7 mb, nf = 2.3 



fm 



leaving the parameters {pDnn and ag^ unchanged. Since the critical density for ip' is 
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Figure 3: (p^) /or J/-?/' as a junction of transverse energy Ef. Data are from for {pi). The 
solid lines are calculated within the threshold model hut with a continuous suppression in the 
time interval [to,''^i]) c./. eq. / fTT]) . 

smaller than that for J /if), more ^'s leak out of the anomalous suppression region, the change 
in suppression due to the increase of is noticeable. Fig. 2 shows the results with a good fit 
to the suppression data. 

The data for {pl)"^' have rather large error bars. The calculated curves show a strong 
dependence on t^i, again with a trend that {pf) at larger values of Et require larger values of t^- 
However, the numerical values tA{Et) for ip' are above those for the J/ip by about 1 — 2/m/c. 
This result is strange, because we expect the ip' to be destroyed more easily and therefore more 
rapidly. We will come back to this point in the conclusions. 

The time structure introduced via eq. (jH]) is certainly oversimplified. Rather than having it 
act at one particular time tA, it is more reasonable to assume that it acts during a time interval. 
Therefore we have also investigated the following form of the absorptive term 



a{s,pt,t; h) 



-e{np{h,s)-n,)e{ti-t)e{t-to). 



ti — to 

Anomalous suppression acts for times t between to and ti (which may be a function of Et). 
Especially to is an interesting quantity, because it gives the starting time for anomalous sup- 
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pression. 

We calculate the suppression S'^^'^{Et) and {pf)^^^'^{Et) as a function of [t^^i] and for ao = 
5. All other parameters are unchanged. The suppression is well described for all intervals 
[to,ti], but the calculations for {p'f)'^^'^(Et) depend strongly on the choice of this interval. Fig. 3 
shows some representative examples: The four windows in Fig. 3 display events with to = 

0. 1.2.3/m/c, respectively, and several values for ti. It is obvious that there is not one curve, 
which describes all the data. Rather we find that the best curve for Et < 60 GeV corresponds 
to [0, 1], for 60 <Et< 80 GeV to [1, 3], for 80 < Et < 100 GeV to [2, 5], and for Et > 100 GeV 
to [3, 5]. The length ti — to of the interval remains approximately constant, while the beginning 
time to increases with Et. We will discuss the significance of these results after we have treated 
leakage within the comover model. 

We proceed to model II, the scenario of comovers: Partons and/or hadrons which move 
with the ip may destroy the ip with a cross section af^. The comover density ricoit) depends 
on time, for which the Bjorken scenario of longitudinal expansion predicts . The absorptive 
term a in eq. (j21) then takes the form 

al.-pl, t; b) = ^^^^!!£M^(!^t„ _ t)e{t - to), (12) 

t Uf 

1. e. absorption by comovers starts at t = tg and ends at ti, when the comover density nc{b, s) ■ 
to/ti has reached a value Uf independent of b and s. The comover approach contains a definite 
time structure for anomalous suppression and we have not changed it. We also account for the 
transverse energy fluctuations [10^, by replacing Uc by ^^*^ nc, and the transverse energy loss[TT] 
induced by the J/ib trigger, by rescaling < i?^ > by a factor which have been shown to 

Tip 

be significant for the explanation of the sharp decrease of S'^^^{Et) at Et > 100 GeV. In the 
choice of parameters we have followed jHH nco(b,s) = 1.5 np{b,s) with the participant 
density from eq. ((Zj), to = 1 fm/c, nj = 1 fm~^, a'^lf = 4.5 mb, afi^^ = 6 mb, a'^l"^ = 1 mb and 
af^ = 3 mb, no additional parameter has been introduced. 

The calculated Et dependence of the suppression shown in Fig. 4 fits the data acceptably 
well like model I. For (pf), we compare the results of two calculations with the data. The 
dashed lines in Fig. 4 show the calculation leaving out leakage. Formally this limit is obtained 
from eq. Q by setting Vt = in the exponent and in Due to the introduction of the Et 
loss which is necessary to recover the J/ip suppression for large values of Et, the case without 
leakage does not fit the data (p^) even in the domain of low Et values. Only when the leakage 
effect is taken into account, the calculation agrees well with the data. We stress: the calculation 
of {pt){Et) in the comover model is a true prediction in the sense that no parameter is adjusted 
above those which are fitted to the suppression S^{Et). We have also calculated the mean time 
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(t^) for comover action by studying the suppression S'^{Et) as a function of time and taking 
the mean of t with the weight dS"^ [Et) / dt and find 

{t^f) = 3.5fm/c 
(t^') = 3.0fm/c, (13) 

which values include the time to, eq. (fT^ . between the end of normal suppression and the 
beginning of comover action. The values eq. (|T3|l are found to be rather independent of Et. 




E, (GeV) 20 40 60 80 100 120 

E, (GeV) 

Figure 4: Nuclear suppression a'^ /a^^ and (pf) for J/ip (above) and ip' (below) as a function 
of transverse energy Ef. Dotted and solid lines are calculated in the comover model without and 
with considering the leakage effect, respectively. 

We summarize: In this letter we have investigated the influence of leakage on the calculation 
of (pf){Et) for J/ip and ip' produced in Pb-Pb collisions at SPS energies. This effect is closely 
related to time structure of anomalous suppression. The evolution of the ip during anomalous 
suppression including leakage is described within a general transport equation. The formalism is 
applied to two models with rather contradictory underling physical assumptions, the threshold 
(QGP) and the comover models with the following results: 

(i) Calculations within the original models, where leakage is left out, do not describe the 
data for {p^)^{Et), the discrepancy being particularly strong at high values of Et. 
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(ii) Including leakage into the comover model, without changing its structure nor its param- 
eters leads to a good agreement with the data for {p^){Et) for J/ip over the full range of 
values Et. 

(iii) The assumption in the threshold model that anomalous suppression acts instantaneously 
at tA = 0, i.e. right after normal suppression is not supported by experiment. Rather 
for central collisions, the data of (p^) are described best, if one assumes that anomalous 
suppression acts at a time = 3 — 4 fm/c after the nuclear overlap. 

(iv) The situation for the ip' is less clear. While the data for suppression can be fitted by 
properly adjusting the parameters, both models underpredict the data for {pf)'^'{Et). 
This is evident for the comover model. Within the threshold model the values required 
to fit the ip' data are larger than those for the J/ip, which seems unreasonable to us. It 
could be that the error bars on the experiments are too small. 

We conclude: In this letter we have investigated how leakage (escape of high pt ip's) when 
introduced into exciting models of anomalous suppression influences the calculated values of 
{p1){Et). As mentioned already above, there could also be other effects which could influence 
{pD- We discuss them briefly: pt dependence of the mechanism responsible for anomalous 
suppression and transverse expansion of the plasma. Both effects can be treated within the 
transport approach eqs. Q,® by introducing an explicit dependence on pt and a modified 
dependence on t into the function a{s,pt,t]b). A first and schematic investigation on the pt 
has been made by Dingh JT^. One sees without calculation that if da/dpt < 0, i.e. high ptip's 
are absorbed less (whatever the mechanism may be), the calculated values of {p'^){Et) increase, 
thus having the same effect as an increase of the time tA eq. (jH)), when anomalous suppression 
acts. As for the transverse expansion of the plasma, its effect is qualitatively clear: It reduces 
the effect of leakage since it makes it harder for the ip's to escape. Although the qualitative 
changes on (pf) of the two effects are clear, their quantitative treatment needs a rather careful 
study of the underlying physics, necessitates to modify existing models, to change their input 
parameters and to introduce new parameters. This is beyond the scope of this letter. 
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